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12  Feb  99 

Successful  mission 

5300m,  6470m,  7690m,  8930m,  1 1850m,  9230m 

13  Feb  99 

Successful  mission 

7898m,  9315m,  10595m,  11580m,  12540m 

20  Feb  99 

Successful  mission 

21  Feb  99 

Successful  mission 

22  Feb  99 

Successful  mission 

23  Feb  99 

Successful  mission 

February  off  the  south  coast  of  Japan  was  chosen  to  provide  measurement  opportunity  in  the  strongest 
winter  jet  stream  winds  that  occur  in  the  northern  hemisphere.  During  the  28  days  of  Februaiy  that  the 
ARA  EGRETT  aircraft  was  at  Yokota  Air  Base  to  the  west  of  Tokyo,  eight  full  measurement  missions 
were  flown  and  one  partial.  The  altitude  range  over  which  these  measurements  were  made  was  12,537 
km  to  5,294  km.  Flight  operations  were  from  Yokota  AB  because  that  was  the  base  that  United  States 
Forces  Japan  (USFJ)  obtained  permission  from  the  Japanese  Foreign  Ministry  for  EGRETT  operations. 
It  was  also  the  only  USAF  base  in  Korea- Japan  region  which  would  guarantee  EGRETT  hanger  space 
during  the  month  of  February.  This  was  essential  to  keep  the  pressure  probes  from  getting  wet  (rain, 
snow,  condensation  etc)  and  freezing  up  at  high  altitudes. 

Yokota  AB  was  a  good  operating  location  related  to  the  airspace  congestion  in  Japan.  To  make 
measurements  at  multiple  levels  on  each  mission,  EGRETT  needed  to  be  in  restricted  airspace. 

Yokota  AB  is  the  location  of  the  5th  Air  Force  to  which  there  was  a  Japanese  Self  Defense  Force 
(JSDF)  liaison  office.  Through  the  collaboration  of  5AF/DOO  and  the  JSDF  liaison  office,  EGRETT 
obtained  permission  to  operate  in  three  restricted  air  space  regions,  two  managed  by  JSDF  and  one 
managed  by  the  US  Navy,  One  JASDF  area  (G-1)  was  at  38N  and  138E;  the  other  area  (K-1)  was  at 
34N  and  137. 5E  (approximately  440  km  separation).  The  US  Navy  area  was  at  34.3N  and  141E. 

Allowing  28  days  to  obtain  about  fifty  hours  of  measurements  was  just  sufficient.  The  weather  was 
very  cooperative.  Only  three  days  were  lost  to  flying  because  of  snow  or  rain.  Seven  flying  days  were 
lost  to  flu  as  all  six  members  of  the  mission  crew  came  down  with  it  at  the  same  time. 

The  turbulence  instrumentation  on  the  aircraft  was  identical  to  the  one  used  during  the  Australian  study 
in  mid- 1998,  but  had  been  improved  again  in  terms  of  response  time  and  reliability.  Only  very  few 
failures  occurred,  mainly  caused  by  insufficient  GPS  satellite  coverage. 

Flight  and  altitude  tracks  of  these  flights  are  shown  in  Appendix  B. 


2.3  Southern  Australia  (14  Jul  1999  - 15  Aug  1999,  Adelaide) 

Eight  research  missions  were  flown  with  a  total  of  36  hours. 


Date 

Comments 

16  Jul  99 

Instrumentation  test  flight 

28  Jul  99 

Instrumentation  test  flight 

29  Jul  99 

Instrumentation  test  flight 

3  Aug  99 

Instrumentation  test  flight 

5  Aug  99 

Calibration  flight 

6  Aug  99 

Successful  mission 

FL290,  FL300,  FL310,  FL320,  FL270,  FL330,  FL390,  FL365 

11  Aug  99 

Successful  mission 

FL310,  FL450,  FL440,  FL360 

12  Aug  99 

Successful  mission 

FL320,  FL305 

Flight  and  altitude  tracks  of  these  flights  are  shown  in  Appendix  B. 


2,4  United  Kingdom  (1  May  2000  -  7  June  2000,  Boscombe  Downs) 

The  campaign  in  the  UK  was  funded  by  a  research  grant  from  the  UK  National  Environmental 
Research  Council  (NERC)  to  the  University  of  Wales  (UWA,  Principal  Investigators:  Drs.  Jim 
Whiteway  and  Geraint  Vaughan/University  of  Wales).  Data  from  this  campaign  was  available  to  the 
current  AFOSR/AOARD  project.  The  campaign  was  conducted  under  the  acronym  AberEgrett,  short 
for  Aberystwyth  -  Egrett  -  Experiment.  The  main  aim  of  AberEgrett  was  to  investigate  atmospheric 
processes  around  the  Tropopause  (around  9- 12km  altitude  at  that  time  of  the  year)  which  influence  the 
transport  of  trace  gases  harmful  to  the  Ozone  Layer  from  low  level  up  into  the  Ozone  Layer.  Such 
processes  include  gravity  waves,  as  well  as  turbulence  which  may  be  caused  by  breaking  waves  in  the 
lee  of  the  Welsh  mountains  or  vertical  wind  shear,  especially  around  the  high  wind  layers  (jet  streams) 
which  occur  frequently  at  that  altitude. 

In  addition  to  the  standard  turbulence  instrumentation  (BAT-Probes),  the  Egrett  carried  sensors  for 
humidity  and  trace  gases.  Amongst  these  was  a  GC  (gas-chromatograph)  developed  at  Imperial  College 
at  Cambridge  University,  a  TDL  (tunable  diode  laser)  developed  by  UKTNPL  (National  Physics 
Laboratory),  a  fast  Ozone  analyser  developed  at  the  German  Aerospace  Agency  (DLR),  a  frost  point 
mirror  also  from  DLR,  and  an  ultra-fast  temperature  (UFT)  sensor  developed  at  the  University  of 
Warsaw  in  Poland. 

Apart  from  the  aircraft  measurements,  there  were  also  measurements  using  radars,  lidars  and  balloons. 
A  total  of  90  hours  were  flown  during  18  flight  missions. 


Date 

Comments 

4  May  00 

Instrumentation  test  flight  to  13km  altitude 

5  May  00 

Chemistry  flight  (aborted) 

10  May  00 

Gravity  wave  flight 

11  May  00 

Gravity  wave  flight 

12  May  00 

Gravity  wave  flight 

13  May  00 

Chemistry  flight 

15  May  00 

Chemistry  flight 

16  May  00 

Chemistry  flight 

17  May  00 

Gravity  wave  flight 

18  May  00 

Chemistry  flight 

19  May  00 

Gravity  wave  flight  (aborted) 

22  May  00 

Gravity  wave  flight  (aborted) 

24  May  00 

Gravity  wave  flight 

25  May  00 

Gravity  wave  flight 

31  May  00 

Gravity  wave  flight 

2  Jun  00 

Chemistry  flight 

o 


5  Jun  00 

Calibration  and  chemistry  flight 

6  Jun  00 

Gravity  wave  flight 

The  flight  patterns  were  designed  with  regard  to  the  scientific  objectives  and  the  meteorological 
conditions.  The  most  basic  mode  was  to  fly  the  Egrett  at  successive  levels  from  below  the  tropopause 
and  up  to  15  km.  These  flights  passed  over  Aberystv^th  in  a  convenient  pattern  (eg.  race-track  or 
figure-eight).  A  number  of  flight  paths  sampled  the  wave  field  above  the  ocean,  making  them  very 
useful  for  identifying  wave  sources  other  than  mountains.  When  possible,  flight  paths  were  designed  to 
observe  the  distribution  of  gravity  wave  activity  and  tracer  filamentation  with  respect  to  the  jet  stream. 
Another  basic  mode  was  to  focus  more  directly  on  the  tropopause:  flying  slightly  above  and  below,  in 
order  to  concentrate  on  stratosphere-troposphere  exchange  by  gravity  wave  breaking  and  turbulent 
mixing.  Measurements  of  ozone  and  water  vapour  indicated  how  the  wave  breaking  and  turbulence  is 
affecting  vertical  transport.  In  this  regard,  upwind  measurements  by  the  Egrett  were  be  indispensable  in 
the  interpretation  of  vertical  profiles  measured  above  Aberystwyth;  and  also  the  converse. 

Apart  from  the  aircraft  measurements,  the  following  observations  were  taken: 

Radar:  The  NERC  MST  radar  located  just  East  of  Aberystwyth  provided  measurements  of  the  three 
dimensional  motion  field  and  turbulence  up  to  an  altitude  of  20  km.  The  basic  technique  is  to  detect  the 
backscatter  of  VHP  (46.6  MHz)  radar  signals  from  small-scale  temperature  and  humidity  fluctuations. 
Wind  velocity  is  derived  from  the  Doppler  shift  and  the  intensity  of  turbulence  is  determined  from  the 
broadening  of  the  radar  backscatter  spectra.  The  MST  radar  is  also  of  utility  for  studies  of  filamentation 
and  layering:  the  radar  backscatter  power  profiles  can  monitor  the  progression  of  long  lived  static 
stability  anomalies  in  the  lower  stratosphere. 

Balloons:  Temperature  and  wind  radiosondes  were  launched  four  times  each  day  from  the  weather 
station  at  Aberporth  (40  km  from  Aberystwyth)  providing  regular  profiles  of  temperature  and 
horizontal  wind.  Radiosondes  launched  at  Aberystwyth  provided  temperature  profiles  (wind  is 
measured  by  the  radar),  and  profiles  of  ozone 

Lidar  measurements:  The  UWA  differential  absorption  lidar  (DIAL)  was  used  for  continuous 
measurements  of  ozone  vertical  profiles  in  the  tropopause  region  with  temporal  and  vertical  resolutions 
of  30  min.  and  500  m.  This  lidar  is  based  on  a  Nd:  YAG  laser  (fourth  harmonic  output  at  266  nm)  'with 
Raman  shifting  in  D2  and  H2  to  generate  the  transmitted  wavelengths  of  289  nm,  299  nm,  and  316  nm. 
The  resolution  of  these  measurements  is  sufficient  to  observe  the  perturbations  induced  by  gravity 
waves  and  also  the  consequences  of  irreversible  mixing  during  gravity  wave  breaking.  This  lidar  also 
indicated  the  presence  of  cirrus  clouds. 


Flight  and  altitude  tracks  of  these  flights  are  shown  in  Appendix  B. 


3  INSTRUMENTATION 
3,1  BAT-Probe: 


Measuring  turbulence  from  aircraft  is  a  technologically  demanding  task.  It  involves,  in  general,  the 
measurements  of  air  temperature,  pressures  and  the  3D-wind  vector  >vith  high  accuracy  and  at  sampling 
rates  of  tens  of  Hertz.  Fast  sensors  for  air  temperature  and  pressures  are  available,  are  light-weight  and 
do  not  require  much  space  and  electrical  power.  Until  recently,  the  main  problem  was  to  measure  the 
3D-wind  vector  accurately. 

Wind  measurements  from  an  aircraft  require  observation  of  the  wind  velocity  relative  to  the  sensors 
and  the  velocity  of  the  sensors  relative  to  the  Earth.  The  3D-wind  vector  is  then  found  as  the  small 
resultant  vector  sum  of  these  two  large  vector  velocities. 

None  of  the  two  large  vector  velocities  can  be  measured  directly,  but  require  the  measurement  of  a 
number  of  other  parameters,  amongst  them,  the  angles  of  attack  and  sideslip  of  the  aircraft,  the 
aircraft’s  heading  and  its  pitch  and  roll  angles,  3D-accelerations  of  the  aircraft,  its  angular  rates  and 
more.  All  measurements  have  to  be  combined  in  a  very  careful  manner,  without  phase  shifts  or  other 
inconsistencies  between  them,  to  derive  the  two  large  vector  quantities.  For  most  applications,  it  is  also 
essential  to  know  the  position  of  the  aircraft  accurately  at  all  times. 

Traditionally,  combinations  of  a  number  of  different  sensor  systems  were  used  to  measure  turbulence, 
amongst  them  Inertial  Navigation  Systems  (INS),  Doppler  systems,  smaller  gyro  systems,  GPS  attitude 
systems,  differential  pressure  port  systems  at  the  nose  section  of  research  aircraft  or  moveable  or  rigid 
vanes  mounted  to  nose  booms.  Most  of  these  systems  involved  one  or  several  very  complex,  heavy, 
power-hungry  and  expensive  components  and  were  usually  only  ‘one-off  dedicated  installations. 

The  BAT-Probe  (it  got  its  name  from  its  shape  -  it  looks  like  a  baseball  bat-  but  it  also  could  mean 
“Best”  Aircraft  Turbulence  Probe)  was  developed  as  a  means  to  overcome  most  of  the  above 
mentioned  complexity  of  installation.  It  is  an  attempt  to  make  the  complex  measurements  basically 
independent  of  the  platform  by  combining  the  latest  sensor  technology  with  powerful  electronics  and 
computing  hard-  and  software  into  a  small  package. 

The  BAT-Probe  combines  (low-cost)  sensors  for  pressure,  air  speed,  angles  of  attack  and  sideslip,  3D- 
acceleration  and  air  temperature  with  differential  GPS  systems  for  high  precision  position,  ground 
speed  and  attitude  (aircraft  pitch  and  roll  angles  and  heading)  measurements  packaged  together  in  a 
small  stand-alone  package.  The  probe  contains  its  own  16-bit  A/D-unit  which  converts  the 
measurements  to  digital  values  directly  at  the  sensor  location  and  transmits  the  data  in  a  serial  data 
stream  to  a  PC-based  real-time  processing  computer. 

Combining  the  sensors  at  one  location  in  close  proximity  of  the  A/D-unit  minimises  possible 
attenuation  of  the  small  signal  voltages.  The  small  size  of  the  probe  with  its  very  short  pressure  lines 
between  the  pressure  ports  and  the  pressure  transducers  also  virtually  eliminates  problems  with  time 
lags  or  resonance  effects  in  pressure  lines. 

The  sensor  package  itself  was  mainly  developed  under  NOAA/ATDD  responsibility,  while  ARA  and 
FIAMS  (the  predecessor  of  ARA)  developed  the  A/D-unit  and  associated  digital  electronics.  Both 
organisations  have  joint  responsibility  for  the  low-  and  high-level  software  development. 

The  probe  weighs  about  4kg  and  requires  about  50  W  of  10  to  30  VDC  power.  Like  a  sonic 
anemometer,  data  is  received  serially  at  50  Hz. 

The  housing  consists  of  a  15-cm  diameter  hemisphere  supported  by  a  truncated  circular  cone.  Its 
carbon-fiber  construction  is  light  but  very  strong.  Further,  carbon-fibers  attenuate  any  potential 
electromagnetic  interference. 

Mechanical  and  electrical  installation  on  any  aircraft  is  simplified  by  placing  electronics  and  sensors 
within  the  housing.  Relative  winds  are  computed  from  the  pressure  distribution  observed  at  nine 
pressure  taps  on  the  sphere.  Two  unique  features  of  the  BAT  probe  are  the  method  of  measuring  static 
pressure  and  atmospheric  temperature. 

Static  pressure  is  difficult  to  measure  on  an  aircraft  because  the  sensed  pressure  is  easily  contaminated 
by  flow  distortion  about  the  static  ports  and  airframe.  The  BAT  probe  mitigates  static  error  by  its 
unique  9-hole  design.  With  the  BAT  probe,  static  pressure  is  obtained  from  a  pneumatic  average  of  the 
four  pressure  taps  placed  41 .8°  aft  of  the  Design  Stagnation  Point  (DSP)  and  45°  off  vertical.  At  zero 
angle  of  attack  (when  the  instantaneous  stagnation  point  is  coincident  with  the  DSP),  the  BAT  probe 


correctly  senses  static  pressure.  However,  like  a  conventional  Pitot  probe,  static  pressure  is  increasingly 
underestimated  with  increasing  angles  of  attack.  With  this  approach,  the  static  error  is  less  with  the 
BAT  design  and  an  analytical  model  exists  to  remove  the  error. 

Solid  state  precision  compensated,  low  pressure  sensors  form  SenSym  (for  example  model  SCXL004) 
or  Data  Instruments  (Model  series  XCX)  are  used.  Such  pressure  sensors  have  a  frequency  response 
approaching  IkHz  with  an  accuracy  of  around  0.2%  of  their  full  scale  reading.  Furthermore,  the  sensors 
are  heated  to  reduce  temperature  related  drift. 

The  actual  sensor  choice  is  optimized  for  the  design  flight  speed.  As  a  result,  sensors  with 
unnecessarily  large  ranges  are  not  used  and  accuracy  within  the  desired  sensing  range  is  high.  For 
example,  a  ±12hPa  sensor  is  usually  used  on  angle  of  attack  and  sideslip  pressure  ports.  This  yields  an 
accuracy  of  ±0.05hPa.  Such  accuracy  is  not  possible  with  more  conventional  aircraft  sensors. 

At  aircraft  flight  speeds,  the  observed  temperature  is  contaminated  by  compressibility  effects. 
Correction  of  the  sensed  temperature  is  functionally  dependent  on  dynamic  and  static  pressure.  The 
BAT  probe  places  a  micro-bead  sensing  element  within  the  DSP  where  both  dynamic  and  static 
pressure  are  known  accurately  and  at  high  frequency.  The  DSP  port  is  designed  to  ventilate  the  micro 
bead  at  lOm/s  with  small  thermal  housing  influence.  The  resulting  housing  sensor  combination  has  a 
0.07s  time  response  that  is  much  faster  than  more  conventional  housings/sensor  combinations. 

For  fast  temperature  measurement,  a  Victory  Engineering  Corp's  0.005"  micro-bead  thermistor  is  used. 
Within  the  BAT  housing,  they  have  a  time  response  of  0.07s  that  is  ten  times  faster  than  a  fast 
Rosemount  sensor  housing  combination.  On  the  other  hand,  the  BAT  probe  temperature  sensor  is  more 
fragile.  Although  the  housing  design  is  great  for  fast  response  and  observing  the  micro-beads  pressure 
environment  for  use  in  temperature  recovery  computations,  it  offers  the  micro-bead  little  protection 
from  foreign  object  damage  (bugs,  rain  drops  etc.).  Regardless,  this  has  not  been  a  serious  problem  as 
micro-beads  are  infrequently  damaged.  To  further  mitigate  this  concern,  the  BAT  probe  electronics 
allow  installation  of  twin  micro-beads. 

The  BAT  probe  relies  strongly  on  DGPS  technology  and  does  not  require  any  gyro-based  systems.  In 
the  past  few  years,  Differential  GPS  (DGPS)  techniques  have  developed  into  a  powerful  measurement 
utility.  Currently  DGPS  measurement  techniques  allow  dynamic  observation  of  position,  velocity  and 
attitude  at  lOHz  and  with  accuracy  of  ±0.5m,  ±2cm/s  and  ±0.05°  respectively.  If  necessary, 
accelerometers  can  be  used  to  extend  these  observations  to  higher  frequencies. 

Position  is  required  to  document  aircraft  location  and  is  not  needed  for  wind  computation.  Regardless, 
simple  DGPS  allows  documentation  of  position  to  ±3m.  Where  necessary,  extension  to  high  accuracy 
sub-meter  techniques  is  not  difficult. 

Using  Doppler  DGPS  techniques,  observation  of  probe  velocity  to  20Hz  and  at  an  accuracy  of  ±2cm/s 
is  simple.  A  main  advantage  of  DGPS  approaches  is  that  one  obtains  the  velocity  of  the  probe  in  Earth 
coordinates.  Again,  accelerometers  within  the  BAT  probe  can  be  used  to  extend  velocity  frequency 
response  to  50Hz. 

Attitude  is  important  since  winds  relative  to  the  BAT  probe  must  be  rotated  to  Earth  coordinates. 
Attitude,  at  lOHz  and  with  an  accuracy  of  ±0.05°,  is  obtained  from  a  Trimble  TANS  VECTOR  GPS 
attitude  receiver.  Although  lOHz  may  be  fast  enough  for  large  aircraft,  it  is  not  for  smaller  aircraft  with 
less  inertia.  For  smaller  craft,  GPS  attitude  output  is  low-passed  then  mixed  with  high-passed  twice 
integrated  differential  acceleration  observations.  The  result  is  attitude,  at  50Hz  and  with  an  accuracy  of 
±0.05°.  Four  small  GPS  patch  antennas  must  be  installed  for  use  of  a  TANS  VECTOR  system. 

The  A/D  system  for  the  BAT-Probe  was  developed  by  FIAMS/ARA  to  provide  high  speed,  high 
resolution,  multi-channel  data  logging.  The  system  consists  of  two  boards  with  an  8  pole  low  pass 
(anti-aliasing)  filter  for  each  of  eight  channels  combined  with  an  8  channel  multiplexer  and  an 
intelligent  16  bit  A/D  board.  The  three  boards  together  comprise  a  signal  conditioned,  16  channel,  16 
bit  A/D  system  with  serial  output.  A  custom  ISA  serial  board  provides  an  intelligent  interface  to  a  PC. 

The  system  is  unique  in  its  combination  of  high  speed,  high  resolution,  high  channel  count,  and  on 
board  intelligence.  Software  filters  and  other  operations  may  be  implemented  on  the  A/D  board, 
relieving  the  PC  processor  from  such  simple,  but  time  consuming  tasks.  Additionally,  multiple  systems 
may  reside  on  one  PC,  allowing  for  very  high  chaimel  count  systems. 

The  hardware  and  firmware  are  configured  to  provide  50Hz  data  output  from  the  16  channels.  Each  of 
the  16  channel  is  first  conditioned  by  the  8  pole  30Hz  low  pass  (anti-aliasing)  filter.  Then  it  is 


interrogated  or  over  sampled  32  times  by  the  16  bit  A/D.  The  resulting  data  is  averaged  and  then 
serially  transmitted  and  timed  tagged  to  ±lms  based  on  the  Isec  pulse  available  from  the  GPS  receiver, 

3.2  Global  Positioning  System 

Two  GPS  receivers  are  installed  in  the  aircraft’s  mission  bay,  a  12-channel  NovAtel  system  for  position 
and  velocity  measurements  using  differential  technology  and  a  Trimble  TANS  Vector  system  for 
aircraft  attitude  (pitch,  roll  and  heading). 


3.3  Data  acquisition  and  real-time  monitoring  system 

The  data  system  installed  in  the  Egrett  consisted  of 

•  a  19”  rack  EMI-proof  Pentium  PC  with  dedicated  EISA  interface  cards  for  communication  with 
the  REM-Boxes  (remote  analogue/digital  converter  units),  MUX-box  (digital  multiplexer  unit  for 
REM-boxes),  the  BAT-Probes  and  the  two  GPSs 

•  a  166MHz  Pentium  CPU,  64Mb  of  memory,  a  2Gb  hard  disk,  an  Ethernet  networking  card 

•  a  TFT  touch  screen  in  the  rear  cockpit 

•  a  small  keyboard  with  a  built-in  trackball 

•  a  built-in  IOMEGA  ZIP-drive  for  data  exchange 

The  PC  runs  under  the  Windows  ’95  operating  system.  The  FIAMS/ARA  real-time  data  acquisition 
program  DAMS  and  the  testing  program  REMTEST  are  installed  and  used  to  collect  the  data  and  to 
test  the  installation.  The  system  connects  directly  to  the  pre-installed  wiring  in  the  aircraft. 


4  DATA  PROCESSING 

Most  of  the  data  processing  was  carried  out  using  the  ARATFIAMS  software  package  RAMF/ARAMF. 

The  original  version  of /MA/F  zur  Auswertmg  Meteorologischer  Forschungsfluge, 

Routines  for  the  Processing  of  Meteorological  Research  Flights)  was  developed  by  J.M.  Hacker  at 
the  Meteorological  Institute  of  the  University  of  Bonn/Germany,  to  process  data  from  research  flights 
with  three  motor-gliders  owned  by  the  DLRs  (German  AeroSpace  Agency)  Institute  for  Atmospheric 
Physics.  Today  it  is  used  at  a  number  of  sites  to  process  data  from  meteorological  research  aircraft  as 
well  as  a  multitude  of  other  more  general  time  series  applications.  The  core  site  for  RAMF  is  ARA. 
ARAMF  is  an  enhanced  version  of  RAMF,  meaning  “Advanced  RAMF”.  More  details  about  ARAMF 
can  be  found  at  website  http://www.AirbomeResearch.com.au. 


5  SOME  DETAILED  RESULTS 

The  strongest  turbulence  of  all  measurement  campaigns  was  found  during  the  flight  on  6  August  1999. 
The  results  shown  relate  mainly  to  this  case.  The  data  obtained  on  the  flights  on  25  and  26  August  1998 
are  also  discussed  in  more  detail.  For  all  other  flights,  a  summary  is  presented  only. 

5.1  Flight  25  August  1998 

5.1.1  Synoptic  situation 

On  the  days  before  25  August  1998,  a  jet  stream  with  its  maximum  located  immediately  above 
Perth/Westem  Australia  was  observed  with  a  maximum  wind  speed  of  approximately  80m/s.  The  jet 
stream  core  moved  east  along  the  south  coast  of  Australia  and  passed  to  the  north  of  Adelaide.  A  low 
pressure  system  formed  south  east  of  Australia  but  did  not  fully  develop.  A  ridge  of  high  pressure  was 
situated  over  the  Great  Australian  Bight. 

The  500hPa  analysis  for  25  August  (Figure  1)  shows  a  low  pressure  trough  over  South  Australia  and  a 
strong  low  pressure  cell  over  Tasmania. 


The  amplitude  of  the  planetary  wave  in  the  polar  front  was  forecast  to  intensify  and  a  new  low  pressure 
system  formed  on  26  August.  The  aviation  weather  forecast  published  on  25  August  at  07:31  UTC 
indicated  a  weak  front  approaching  Adelaide  at  04:00  UTC. 

Broken  altocumulus  and  altostratus  were  expected  between  FLIOO  (~3km)  and  FL200  (~6km)  with 
broken  stratus  between  1,500ft  ('-0.5km)  and  2,500ft  (0.8km)  with  rain  east  of  the  front  and  with 
showers  west  of  the  front.  To  the  west  of  the  front  scattered  cumulus  was  expected  between  2,500ft 
(0.8km)  and  10,000ft  (--3100). 


Figure  3:  Longitudinal  cross-section  of  the  temperature  field  at  SOOhPa  (left)  and  200hPa  (right). 
Data  shown  is  1-day  mean  for  25  August  1998.  (from  NOAA  website  http://www.noaa.gov) _ 


At  the  time  of  the  flight,  the  edge  of  the  jet  stream  was  overhead  Adelaide.  As  is  shown  in  the  Adelaide 
Airport  sounding  in  Figure  2,  wind  speed  increased  from  just  15m/s  at  5,500m  to  75m/s  at  9,000m  with 
a  moderate  decrease  above.  The  wind  profile  was  linked  to  the  large  temperature  gradient  of 


approximately  2K/100km  between  cold  air  to  the  south,  east  of  Adelaide,  and  a  warm  continental  air 
mass  north  west  of  Adelaide.  Figure  3  shows  a  longitudinal  cross-section  of  the  temperature  field.  Note 
the  reversal  of  the  temperature  gradient  overhead  Adelaide. 


5.1.2  Measurements 

The  flight  on  25  August  1998  commenced  at  18:1  ILT  and  ended  at  20:51LT.  Turbulence  was  found  at 
around  7,000m,  10,000m  and  1 1,800m  altitude.  The  mean  wind  direction  was  215°true.  Seven 
horizontal  legs  were  flown  between  FL360  (~1 1km)  and  FL460  ("-14km),  see  Figure  4. 

After  climbing  up  to  FL360  (~1 1km),  the  Egrett  headed  into  the  wind  and  flew  legs  at  FL360,  FL380, 
FL400,  FL420,  FL460,  FL400,  FL360,  ie.  11.0km,  11.6km,  12.2km,  12.8km,  14.0km,  12.2km,  11km 
altitude). 


The  flight  legs  relative  to  the  temperature  and  wind  speed  field  are  shown  in  Figure  5  which  was 
constructed  from  routine  soundings  at  Woomera,  Ceduna,  Adelaide  and  Mt.  Gambier  (note  the  reversed 
abscissa  with  East  left  and  West  right).  The  jet  core  was  located  at  12km  altitude  just  below  the 
tropopause  to  the  north  west  of  Adelaide  with  wind  speeds  of  approximately  50m/s.  The  largest 
temperature  gradient  can  be  found  just  south  of  Adelaide. 

Vertical  profiles  of  temperature,  dewpoint,  wind  speed,  wind  direction  and  potential  temperature  as 
measured  during  the  descent  of  the  Egrett  are  shown  in  Figure  6.  Turbulence  was  computed  as  the 
RMS  of  the  3D-accelerations.  The  jet  stream  was  approximately  2,500m  thick  with  a  maximum  shear 
of  17kt  per  1,000ft  (2.8m/s  per  100m).  A  significant  change  of  wind  direction  together  with  strong 
wind  shear  lead  to  a  thin  layer  of  turbulent  air  at  7km  altitude.  The  temperature  in  this  layer  was 
constant  with  an  inversion  at  the  top  of  the  turbulent  layer  of  the  order  of +1K.  More  turbulence  was 
found  at  1 1km  and  14km  altitude. 


20.00 


Figure  5:  The  large  diagram 
shows  vertical  cross  sections  of  air 
temperature  in  (solid  lines)  and 
wind  speed  in  m/s  (broken  lines) 
as  constructed  from  radiosonde 
soundings  at  Woomera,  Ceduna, 
Adelaide  and  Mt  Gambier  at 
OOUTC  on  25  August  1998.  Also 
shown  is  the  location  of  Adelaide 
(green)  and  the  altitudes  of  the 
horizontal  flight  legs  (blue).  Note 
the  reversed  abscissa  with  East  left 
and  West  right. 
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The  small  diagram  shows  the 
location  of  the  cross-section  (red), 
the  flight  track  (green)  and 
coastline,  as  well  as  the  1-day 
mean  of  the  SOOhPa  wind field. 
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Figure  6:  Vertical  profiles  of  atmospheric  parameters  as  measured  by  the  Egrett  during  the  descent  of 
the  flight  on  25  August  1998.  Shown  are:  Turbulence  (see  text)  in  arbitrary  units  (yellow);  wind  speed 
in  m/s  (grey);  wind  direction  in  °true  (scaled  by  a  factor  of  10,  magenta);  potential  temperature  in  °C 
(brown);  and  air  temperature  in  (light  blue). 
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5.2  Flight  26  August  1998 
5.2.1  Synoptic  Situation 

As  the  jet  core  of  25  August  1998  moved  east,  a  pool  of  cold  air  remained  near  Adelaide.  Another  jet 
stream  approached  from  Western  Australia  but  was  too  far  away  to  influence  the  upper  airflow  in  South 
Australia.  The  approaching  front  weakened  and  a  new  low  pressure  cell  formed  south  west  of  Perth. 

The  high  pressure  cell  moved  south  to  Tasmania.  The  500hPa  chart  (not  shown)  for  26  August  1998 
showed  a  high  pressure  ridge  over  Perth  and  a  trough  over  New  South  Wales.  The  wind  speed  in  the  jet 
stream  was  expected  to  increase  to  approximately  45m/s  in  the  early  hours  of  26  August  1 998  with  a 
maximum  increase  of  wind  speed  of  about  20m/s  between  5.6km  and  7.3km  ('^1.2m/s  per  100m).  A 
wide  horizontal  band  of  more  or  less  constant  wind  speed  was  then  forcast  for  the  next  six  hours.  The 
radiosonde  data  indicate  33m/s  wind  speed  at  the  400hPa  level,  the  bottom  of  the  jet  stream,  and  34m/s 
at  the  150hPa  level,  the  top  of  the  jet  stream. 

A  front  passed  Adelaide  the  day  before  and  brought  a  few  showers  and  the  temperature  decreased  by 
about  4K.  As  this  system  passed,  the  strong  temperature  gradient  in  the  area  above  Adelaide  weakened 
in  the  upper  troposphere.  The  cold  air  moved  east  and  warmed,  as  can  be  seen  in  Figure  7  when 
comparing  it  to  Figure  3, 


Figure  7:  Longitudinal  cross-section  of  the  temperature  field  at  SOOhPa  (left)  and  200hPa  (right). 
Data  shown  is  1-day  mean  for  26  August  1998.  (from  NOAA  website  http://www.noaa.gov) 


Figure  8  shows  the  250hPa  winds.  Note  the  confluence  zone  near  Adelaide. 


5.2.2  Measurements 

The  flight  on  26  August  1998  commenced  at  17:42LT  and  ended  at  21:37LT.  Nine  horizontal  legs  were 
flown  between  FL210  ('^6.4km)  and  FL435  ('^1 3.3km),  SQQFigure  9. 

After  climbing  up  to  FL360  (~1 1km),  the  Egrett  headed  into  the  wind  and  flew  several  legs  at  FL210, 
FL360,  FL380,  FL400,  FL420  and  FL435,  ie.  6.4km,  1 1.0km,  1 1.6km,  12.2km,  12.8km  and  13.3km 
altitude). 


Figure  9:  Horizontal  and  vertical  flight  path  for  flight  on  26  August  1998 


The  flight  legs  relative  to  the  temperature  and  wind  speed  field  are  shown  in  Figure  10  which  was 
constructed  from  routine  soundings  at  Esperance,  Forrest,  Ceduna,  Adelaide  and  Wagga  Wagga  (note 
the  reversed  abscissa  with  East  left  and  West  right). 
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Figure  10:  The  large  diagram  shows  vertical  cross  sections  of  air  temperature  in  °C  (solid  lines)  and 
wind  speed  in  m/s  (broken  lines)  as  constructed  from  radiosonde  soundings  at  Esperance,  Forrest, , 
Ceduna,  Adelaide  and  fVagga  Wagga  at  OOUTC  on  25  August  1998.  Also  shown  is  the  location  of 
Adelaide  (green)  and  the  altitudes  of  the  horizontal  flight  legs  (blue).  Note  the  reversed  abscissa  with 
East  left  and  West  right.  The  small  diagram  shows  the  location  of  the  cross-section  (red),  the  flight 
track  (green)  and  coastline,  as  well  as  the  1-day  mean  of  the  SOOhPa  wind  field 
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Figure  11:  Vertical  sounding  with  the  Egrett  aircraft  on  26  August  1998 


For  26  August  1998,  there  was  no  data  available  from  a  routine  sounding  at  Adelaide  airport. 

Figure  1 1  shows  the  Egrett  data  for  the  descent  at  21 :37LT.  The  thick  layer  of  wind  speed  greater  then 
30m/s  as  forecast  by  the  Bureau  of  Meteorology  was  not  confirmed  by  the  aircraft  measurements. 

The  jet  stream  had  a  narrow  maximum  at  1 1km  altitude  with  wind  speeds  up  to  34m/s  and  a  turbulent 
layer  about  250m  below  this  maximum.  The  wind  direction  stayed  constant  throughout  the  upper 
troposphere  at  around  170°  to  180°true. 


5.3  Flight  6  August  1999 
5.3.1  Synoptic  Situation 


The  situation  in  the  preceding  days  indicated  a  strong  jet  stream  extending  from  Western  Australia 
towards  the  Great  Australian  Bight,  heading  east.  The  24  hour  forecast  of  NCEP  indicated  wind  speeds 
of  between  90  and  lOOm/s  for  8  August.  The  wind  speed  was  forecast  to  decrease  during  the  day  to 
between  80  and  90m/s  whilst  it  passed  the  Great  Australian  Bight.  Adelaide  was  at  the  diffluent  exit 
area  of  the  jet  with  decreased  wind  speeds  of  between  70  and  80m/s.The  grid  point  forecast  from  the 
Bureau  of  Meteorology  indicated  a  similar  situation.  It  positioned  the  jet  maximum  at  about  500km 
south  of  the  Great  Australian  Bight  with  a  maximum  speed  of  155kt  (75m/s)  at  the  200hPa  level  and 
with  a  cyclonic  curvature  to  the  north  in  the  vicinity  of  Eyre  Peninsula. 

As  observed  in  the  preceding  weeks,  the  jet  was  expected  to  pass  Adelaide  to  the  north  and  out  of  range 
of  the  aircraft.  An  inspection  of  the  mean  sea  level  chart  from  the  Bureau  of  Meteorology  indicated  a 
partially  occluded  but  weak  cold  front  system  centered  south  of  Perth,  with  a  large  meridional  extent 
from  latitude  45°S  to  23°S.  Adelaide  remained  cloud  free  with  some  advective  cumulus  clouds  up  to 
1,000m  and  excellent  visibility.  Wave-like  cirrus  clouds  were  observed  in  narrow  ranges.  The 
tropopause  height,  as  given  by  the  Bureau  of  Meteorology,  was  1 19hPa  (15,200m)  with  a  temperature 
of-63.rC. 

Figure  12  shows  the  two  closest  aerological  diagrams  of  the  flight  track  in  terms  of  time.  The  data  was 
taken  from  the  radiosonde  soundings  at  Adelaide  Airport.  The  green  and  cyan  lines  are  the  temperature 
and  dewpoint  temperature  respectively,  with  the  scale  on  the  top  axis,  whilst  the  blue  line  is  the  wind 
speed.  The  wind  direction  is  indicated  by  the  magenta  line  and  scaled  by  division  of  10  around  north  (0 
on  bottom  axis).  The  radiosonde  released  at  23:00UTC  (08:30LT)  indicated  very  dry  conditions  above 
the  boundary  layer  up  to  1,000m  due  to  a  nocturnal  inversion.  The  temperature  increase  in  the  inversion 
was  approximately  2.3°C  in  a  layer  approximately  140m  thick.  A  second  inversion  at  9,500m  and  a 
strong  change  of  wind  direction  between  7,500  and  8,500m  was  also  evident. 

The  first  significant  change  was  at  the  height  of  1,000m,  with  an  inversion  on  top  of  the  boundary  layer 
due  to  radiative  cooling  during  the  night.  A  few  clouds  were  observed  just  below  the  inversion.  The 
temperature  then  followed  the  temperature  gradient  of  the  standard  atmosphere  with  dry  to  very  dry 
conditions  throughout  the  middle  troposphere.  The  average  wind  speed  up  to  7,000m  was  lOm/s 
decreasing  to  5m/s  at  8,000m  combined  with  a  rapid  change  of  direction  from  40°true  to  240®true 
which  was  the  direction  of  the  jet  stream. 

The  boundary  layer  height  increased  to  1,500  m  during  the  day  due  to  heating  of  the  Earth‘s  surface. 
Very  dry  conditions  were  found  above.  A  layer  of  much  moister  air  was  advected  from  the  west.  This 
layer  was  about  2,000m  thick  topped  by  an  inversion  of  aproximately  +2K.  The  overlying  airmass  was 
very  dry.  The  wind  speed  increased  during  the  day  because  of  the  jet  stream  maximum  which  passed 
Adelaide  to  the  north  during  the  day. 
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The  mean  sea  level  analysis  at  09:3LT  (Figure  13)  indicated  that  the  cold  front  was  particularly  weak 
with  the  center  of  the  low  pressure  system  to  the  south  of  Albany  moving  east  at  approximately  15m/s. 
The  jet  axis  is  shown  as  a  blue  line.  The  maximum  wind  speed  was  85m/s  directly  above  the  cold  front. 
The  system  was  becoming  occluded  on  the  next  day  with  the  wind  speed  decreasing  to  lOm/s  (see 
Figure  14).  The  24  hour  isallobaric  analysis  (Figure  15)  supported  the  indications  about  the  weakness 
of  the  cold  front  with  a  maximum  pressure  change  of  -2.7hPa  per  day.  The  next  fast  moving  low 
pressure  system  was  approaching  from  the  south  so  that  the  high  pressure  area  was  forced  to  move 
north  east. 


Figure  16:  Max  wind  and  Tropopause  analysis  6  August  1999  OOUTC 


The  Bureau  of  Meteorology* s  maximum  wind  and  tropopause  analysis  from  6  August  1999  shows  the 
wave-like  axis  of  the  jet  stream  and  upper  atmospheric  temperature  conditions  {Figure  16).  Further 
information  about  the  synoptic  conditions  can  be  taken  from  Figure  17  and  Figure  18. 
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Figure  1 7:  Wind  and  temperature  chart  of  the  Australian  region  valid  for  6  August  1999. 
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Figure  18:  Grid  point  forecast  for  the  Australian  region  valid  for  6  August  1999  06UTC. 


5,3.2  Measurements 


The  flight  on  6  August  was  one  of  the  most  successful  flights  in  this  project.  Very  strong  turbulence 
was  found  in  a  thick  layer  between  8,000m  and  9,200m  with  total  acceleration  sometimes  exceeding 
1.5g. 

The  flight  commenced  at  13:39LT  and  ended  at  18:44LT  {Figure  19).  After  climbing  to  29,000ft 
(8.8km),  the  aircraft  was  heading  straight  into  the  wind,  flying  quasi  horizontal  levels  at  FL290,  FL300, 
FL310,  FL320,  FL365  and  FL390,  ie.  8.8km,  9.1km,  9.5to,  9.8km,  11.1km  and  11,9km,  Each  leg  was 
approximately  30min.  In  between  there  were  two  legs  flown  from  FL320  (11.1km)  down  to  FL270 
(8,2km)  and  back  up  to  FL330  (10.0km). 


no 


alt  [m] 


flight  down  990806 


-65  -60  -55  -50  -45  -40  -35  -30  -25  -20  -15  -10  -5  0  5  10  15  20  25 


~ 

rH 

n 

h 

— 

+- 

r-H 

M 

1 

rh 

T- 

r+i 

■l-l 

n 

r+i 

r 

rh 

:5 

0 

L_ 

1 

7 

L 

L 

t 

/ 

4J3 

•  Oi 

r 

1 

I 

1 

■ 

■ 

■ 

■1 

1 

fi 

■ 

■ 

■ 

■ 

■ 

■ 

i! 

■ 

■ 

i 

1 

■ 

1 

I 

1 

■ 

■ 

■ 

■ 

■ 

B 

SI 

■ 

■ 

\ 

■ 

■ 

■ 

E 

■ 

II 

■ 

■ 

■ 

g 

S 

■ 

1 

■ 

■ 

1 

■ 

■ 

■I 

i 

H 

9 

_ 

— 

7 

" 

P! 

1 

i 

■ 

■ 

■1 

1 

_ 

■ 

■ 

IS 

i 

■ 

■ 

■ 

1 

■ 

Si 

-JSt 

LiJ  .'v 

■ 

B 

■ 

1 

■ 

■ 

5 

■ 

1 

1 

1 

^1 

y 

L. 

■ — 

/ 

\ 

> 

7 

s. 

i 

tcs:: 

7*^ 

I 

r-'- 

_ 

c 

V '  '‘H 

L 

-.*1. 

y 

t: 

/ 

c 

-4 

1 

7 

IL 

Z 

1 

r 

i- 

t 

~7 

t 

....'^ ... 

' 

_ _ _ 

7^ 

z 

I 

1 

1 

i^r 

1 

'-V 

\ 

1 

7 

J 

_i 

3- 

1 

r 

7 

ia 

44 

_ _ 

j 

;l_ 

-15  -10  -5 


10  15  20  25  30  35  40  45  50  55  60  65  70  75  80  85  90  95 


-  12500 

^2000 

J1500 

-  11000 

-  10500 

-  10000 

-  9500 

-  9000 

-  8500 

-  8000 

-  7500 

. turbulence 

-  7000 

- windspeed  [m/s] 

-  6500 

- winddir[deg  true] 

-  6000 

- potential  temp 

•  5500 

temp 

-  5000 

—  dewpoint 

-  4500 

-  4000 

-  3500 

-  3000 

-  2500 

-  2000 

*  1500 

-  1000 

-  500 

Figure  21:  As  Figure  1 1,  but  for  the  flight  on  6August  1999. 
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APPENDIX  A:  Conference  Paper 


CLEAR  AIR  TURBULENCE  AND 
REFRACTIVE  TURBULENCE  IN  UPPER 
TROPOSPHERE  AND  LOWER 
STRATOSPHERE 


Owen  R.  Cote*,  Air  Force  Research  Laboratory,  Space 
Vehicles  Directorate,  Hanscom  AFB, 

Jorg  M.Hacker,  Airborne  Research  Australia,  Flinders 
University,  5016  Salisbury 

South,  Australia,  Timothy  L.  Crawford, 
NOAA/ERL/ARL,  Idaho  Falls,  and 

Ronald  J.  Dobosy,  NOAA/ERL/ARL,  Oak  Ridge. 


1 .  Introduction. 

Refractive  Turbulence,  the  local  fluctuations  of 
the  propagation  speed  of  electromagnetic 
waves  relative  to  speed  in  a  vacuum,  is 
associated  with  atmospheric  density 
fluctuations  and  is  a  main  contributor  to  fast 
scintillation  of  laser  communication,  wide-area 
microwave  and  infra-red  surveillance,  and 
directed  energy  weapon  systems.  Slow  fading 
and  the  phenomena  of  radar  or  radio  “holes” 
defined  by  periods  of  anomalous  low  signal 
strength  are  associated  with  strong  and 
persistent  mean  vertical  gradients  of  refraction 
that  occur  with  humidity  layering  and 
temperature  inversion  layers.  The  goal  of 
direct  airborne  measurement  of  refractive 
turbulence  is  to  understand  and  hopefully 
predict  the  evolution  of  both  mean  and 
turbulent  atmospheric  structure  that  causes 
adverse  performance  of  systems  such  as  an 
airborne  early  warning  radar  or  an  airborne 
laser  communication  system.  In  the  upper 
atmosphere  the  refractive  turbulence  parameter 
of  primary  concern  is  the  temperature  structure 
constant,  Ct^ 

Because  the  US  Air  Force  concern  is  with 
global  operation  of  such  systems  the  need  is  to 
know  in  what  geographic  regions,  in  what 
season,  and  in  what  part  of  a  large-scale 
meteorological  pattern,  their  operations  could 
be  most  severely  degraded  by  refractive 
turbulence.  To  answer  the  question  on  just  how 
severe  refractive  turbulence  can  be  a 
collaborative  research  program  to  directly 
measure  refractive  turbulence  above  and  below 
the  core  of  the  winter  sub-tropical  jet  stream 
between  the  altitudes  of  8  and  15  km  was 
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developed.  The  GROB  520T  EGRETT, 
owned  and  operated  by  Airborne  Research 
Australia  (ARA)  was  chosen  to  be  the  airborne 
platform  with  the  NOAA/FRD  BAT  probe 
chosen  as  a  turbulence  sensor.  The  BAT  probe 
is  a  13  cm  diameter  pressure  hemisphere  with 
nine  pressure  ports.  A  micro-bead  thermistor 
for  temperature  measurement  is  mounted  in  the 
central  pressure  port.  The  BAT  probe  is 
designed  to  be  flown  on  aircraft  without 
inertial  navigation  systems.  Each  probe 
contains  a  GPS  receiver  and  three  component 
accelerometers.  For  a  more  complete 
description  of  the  BAT  probe  see  Crawford 
and  Dobosy  (1992).  A  more  recent  description 
is  available  on  the  web  site: 
http://www.noaa.inel.gov/frd/Capabilities/BAT 
/. 

Three  of  the  BAT  probes  were  mounted  on 
EGRETT,  one  under  each  wing  and  one  at  the 
top  of  the  tail.  The  use  of  three  probes  was  a 
response  to  several  problems  relative  to  high 
altitude  turbulence  measurements.  The  first  is 
the  averaging  time  problem  discussed  by 
Lumley  and  Panofsky  (1964),  and  Wyngaard 
(1973).  The  second  problem  is  flow  distortion. 
It  is  hoped  that  the  multiple  probes  will  help  in 
the  identification  of  flow  distortion  affects  on 
spectra.  The  final  problem  is  to  obtain 
measurements  of  vertical  gradients  of  mean 
horizontal  wind  and  mean  potential 
temperature  with  comparable  accuracy  to  the 
turbulence  estimates.  Both  measurements  are 
used  to  evaluate  budgets  of  kinetic  energy, 
temperature  variance,  and  heat  flux. 


2.  Where  and  When  to  Measure  for 
Turbulence. 

Because  funding  was  limited  it  was  important 
that  any  measurement  campaign  should 
provide  information  relevant  to  the  range  of 
maximum  temperature  structure  parameter, 

Ct^,  that  occurs  above  9  km.  The 
determination  of  where  to  measure,  in  what 
season,  and  at  what  altitudes  to  measure  was 
greatly  helped 

by  use  of  climate  data  available  on  the  NOAA 
web  site  (http://www.cdc.noaa.gov/HistData/ ). 
Using  this  data  base  it  was  concluded  that  in 
the  winter  subtropical  jet  stream,  which  is 
found  near  30  degrees  latitude  in  both 
hemispheres,  some  of  the  strongest  winds  are 
found  near  the  south  coast  of  Japan  in  the 
northern  hemisphere  and  in  the  Perth  to 
Adelaide  corridor  in  the  southern  hemisphere. 
Associated  with  the  jet 


core  is  a  change  in  tropopause  heights.  On  the 
pole- ward  edge  it  was  about  9  to  10  km  and 
300  mb  pressure.  On  the  equator-ward  edge  it 
was  varied  from  15  to  16  km  and  about  120 
mb  pressure.  The 

daily  plot  of  the  spatial  distribution  of  the 
tropopause  pressure  jump  clearly  identifies  a 
narrow  ribbon  four  or  five  degrees  of  latitude 
wide  which  migrates  toward  the  pole  in  spring, 
and  retreats  toward  the  equator  in  the  fall 
season.  Some  estimate  of  relative  turbulence 
intensity  can  be  obtained  from  the  temperature 
plots  for  the  standard  pressure  levels.  The 
magnitude  of  horizontal  temperature  gradients 
on  a  constant  pressure  surface  indicates  the 
strength  of  the  vertical  gradient  of  the 
geostrophic  wind.  A  narrow  ribbon  of  the 
maximum  temperature  gradient  at  each 
pressure  coincided  with  the  location  jet  core 
and  the  tropopause  jump  with  the  greatest 
geostrophic  shear  found  above  and  below  the 
jet  core  at  the  150  mb  and  300  mb  pressure 
levels.  Examining  several  years  of  this  led  to 
the  choice  of  February  in  Japan  and  August  in 
Australia. 

Having  made  the  turbulence  measurements  in 
these  months  in  1999,  the  NOAA  database  is 
still  useful  in  estimating  the  order  of 
magnitude  of  the  horizontal  shear  and 
stratification  production  terms  in  the  energy, 
flux,  and  temperature  variance  budgets.  The 
upper  troposphere  and  lower  stratosphere  are 
not  homogeneous.  The  measurement  flights 
where  accordingly  in  the  direction  in  which  the 
atmosphere  is  most  homogeneous,  either 
upwind  or  downwind.  The  length  of  a  constant 
altitude  measurement  was  1200  seconds, 
which  is  about  100  kilometers  at  EGRETT 
airspeed.  In  Japan  the  smallest  separation  of 
measurement  levels  was  660  meters  and  in 
Australia  it  was  330  meters. 


3,  Measurements 

In  Table  1  are  summarized  the  structure 
function  parameters  for  temperature  for 
temperature  and  the  three  velocity  components 
for  most  of  the  measurements  made  in  Japan  in 
February  1999.  Also  listed  in  the  table  are  the 
ratios  of  pairs  of  structure  parameters.  C^/  Ct^ 
is  equal  to  the  ratio  ai  e/Pi  %,  the  product  of 
the  dissipation  ratio  with  the  Kolmogorov 
spectral  constant  ratio.  It  shows  great 
variability  in  a  very  reduced  value.  The 
estimates  for  turbulent  dissipation  (e)  are 
obtained  from  the  relation:  =  2  The 

column  labeled  “dyn  heat”  is  the  ratio  of 
measured  to  true  Gp^  if  corrections  for 


dynamic  heating  are  not  made  to  the 
temperature  time  series.  The  ratio  of  the 
transverse  velocity  structure  parameters  to  the 
longitudinal  velocity  structure  parameter 
should  be  4/3  if  the  inertial  sub-range  is  locally 
isotropic.  On  average  Cy^  is  and  Cw^  is  not. 
This  same  result  was  obtained  for  the 
measurements  in  Australia  except  for  6  August 
1999  which  are  also  shown  in  Table  1.  The 
measurement  duration  at  level  9650  m  was 
3000  seconds  in  contrast  to  the  other  levels 
that  day  which  all  were  for  1200  seconds.  The 
time  series  for  measurement  flight  level  was 
dived  into  two  sections.  Section  A  is  the 
eastern  half  and  B  is  the  western  half.  Level 
9560B  had  the  largest  structure  parameters 
EGRETT  has  measured  to  date. 
Correspondingly  the  vertical  accelerometer 
data  for  section  B  was  intermittently  ±  1.5  g. 
On  the  following  two  days,  7and  8  August,  it 
rained  and  EGRETT  could  not  make 
measurements  because  water  in  the  pressure 
ports  would  freeze  at  altitude.  It  is  of  interest 
to  look  at  the  NOAA  CDC  web  site  (see 
above)  for  these  two  days.  The  presence  of  a 
low-  pressure  trough  (with  closed  temperature 
contours)  contiguous  to  and  pole-ward  of  the 
jet  is  the  structure  we  identified  as  most 
favorable  to  extreme  turbulence  and  we  sought 
to  investigate  whenever  possible. 

In  figure  1  are  vertical  profiles  for  Gp^ ,  one 
from  Japan  and  one  from  Australia  that  have 
the  peak  values  for  each  campaign.  With  the 
660  meter  separation  between  measurement 
levels  in  Japan  it  is  possible  that  the  peak  value 
was  not  sampled.  To  capture  the  structure  of 
these  narrow  layers  of  extreme  turbulence 
future  measurements  should  probably  be 
between  250  and  300  meter  separation  in  the 
vertical  when  extreme  turbulence  is 
encountered.  Figure  2  contains  the  vertical 
profiles  of  the  three  components  of  velocity 
variance,  temperature  variance,  and  the  vertical 
heat  flux  magnitude  (note  log  plot  of 
amplitude).  Figure  3  contains  the  three 
Reynolds  stress  components  in  a  linear  plot 
that  sets  off  the  turbulent  layer  rather  well.  The 
curvature  of  the  <uw>  profile  in  the  turbulent 
layer  acts  as  a  negative  feedback  and 
diminishes  a  positive  mean  vertical  gradient  of 
the  mean  wind  that  produces  negative  <uw>. 

At  low  spatial  wave-numbers,  velocity  and 
temperature  spectra  follow  a  ki'^  power  law 
and  preliminary  analysis  with  the  scaling 
suggested  by  Lumley  (1964)  and  Weinstock 
(1978)  indicates  a  non-constant  spectral 
constant  in  agreement  with  the  aircraft 
turbulence  measurements  of  Lilly  and  Lester 


(1974).  This  is  an  interesting  preliminary  result 
but  needs  more  confirmation. 

In  Table  two  are  show  various  turbulent 
parameter  ratios  and  scale  calculations.  The 
column  <u0>/<w0>  is  thehorizontal  to 
vertical  heat  flux  ratio.  The  column  <0^>/w^> 
is  the  ratio  of  temperature  variance  to  vertical 
velocity  variance.  The  Ri  <w0>  column  is  the 
vertical  heat  flux  Richardson  number 
equivalent  to  i  The  column  labeled 
<0VCt^  is  the  integral  scale  estimate  obtained 
from 

(<0^/Ct^)^^^  .  The  next  column  is  the  integral 
scale  obtained  from  ((u^+  v^  +w^)/2)^^^/e  =qVe. 
The  last  column  labeled  Integral  Scale  Ratio  is 
the  ratio  of  these  two  integral  scales. 

4.  Scales  and  Budgets 
TEMPERATURE  VARIANCE  BUDGET 

3  (0^  )/3 1  =  -  advection  -  (w0)  d  0/3  z  ( 1  + 
Horizontal  Prod/Vertical  stratification  Prod  )  + 
-  3  (w0^)/3  z  (1+  Horizontal  TT/Vertical  TT) 

where  %  (x  =  k(  (  3fr3xi  f)  )  is  dissipation 
rate  of  temperature  variance  and  TT  stands  for 
turbulent  transport  by  third  moments.  The 
temperature  structure  parameter  is  obtained 
from  the  inertial  sub-range  temperature 
spectrum  in  the  form 

OtOci)  =  1/4  Cr^  k,-“ 

Under  Kolmogorov  scaling  Ct^  =  4  pi  %  /e‘^^ 
and  p]  is  the  one  dimensional  Kolmogorov 
spectral  constant.  Under  stationarity, 
horizontal  homogeneity,  and  universal 
equilibrium,  production  and  dissipation  are 
equal  so  that: 

-(w0)30^z  =  %  =  K  <  (30^Xi)^ )  ; 

To  derive  divide  by  4pK  The 
resulting  identity  evolves  as  follows: 

4Pi  <w0)  30«z  =  4Pi  Kh/e^^  (3Q^z)^  - 
Ct^  =  4Pi  <(30^Xi  )^>  =  4Pi  (dOdzf 


mean  gradient  times  a  ratio  equivalent  to  a 
Reynolds  number:  KH/k  (30^z)2  =( 
(3fr3xi  )2).  The  basic  problem  in  modeling 
CT2  in  stable  stratification  is  to  determine  the 
composition  of  L4/3  or  KH  /el/3  which  has 
the  dimension  of  (length)4/3.  There  are  two 
possibilities  suggested  by  the  EGRETT 
temperature  spectra.  L  is  equal  to  the 
temperature  integral  scale,  £ .  The  second 
possibility  for  is  a  combination  of  length 
scales:  =  £tI£  where  t  is  the  Ellison 

scale  defined  by  (0^)  =  I  ^  (3Q^z)^ .  This 
scale  arises  naturally  in  the  vertical  heat  flux 
conservation  equation. 

An  alternative  form  for  modeling  Ct^  is 
obtained  by  replacing  the  vertical  gradient  of 
potential  temperature  by  0o/g  not  replacing 
the  vertical  heat  flux  by  an  eddy  diffusion 
coefficient  as  follows: 


-4p,  /e‘'^  (we>  ae^z  =  -4p,  g/e^  (we)  /e  n" 
(0o/g  f  8^"  =  =  -2p,  g/00  (we)  /e  N'  (0o/g  f 
Cu'  =  Ct^ 


VERTICAL  HEAT  FLUX  BUDGET 

3(we)^t  =  -  advection  -  (w^)  30^z  (1  + 
Horizontal  -  Production  )  +  g/0  . 

Vertical  -  Production 

l/po(03p/3z)  +  2^cos(p  (u0)  -  3(w^0V3z  (1  + 
HorTransport  )  +  (w0)(3U/3x  +  3V/3y) 
VertTransport 


In  the  stable  boundary  layer,  Wyngaard  et.  al. 
(1971)  showed  the  vertical  heat  flux  budget  is 
a  balance  between  a  gain  from  stratification 
production  (-  ( w^ )  303z)  and  loss  from 
buoyancy  (+  g/0  (0^))  and  the  fluctuating 
pressure  term  (-  l/po(03p/3z).  The  ratio  of 
buoyancy  to  stratification  production  is 
equivalent  to  £  jl  £  where  (w^)  / 

=  £  ^ .  For  measurements  at  the  Boulder 
Tower  under  conditions  of  stable  stratification, 
Hunt  et.  al.  (1985)  found  this  ratio  was  0.64 
for  all  cases.  In  the  6  Aug  clear  air  turbulence 
event  this  ratio  was  between  1 .4  and  2.4  for 
the  strong  turbulence  levels.  It  is  not  clear  at 
this  time  what  maintains  the  negative  vertical 
heat  flux. 


In  the  above  identity  we  have  replaced  mean- 
square  fine  scale  gradients  by  the  square  of  the 


Conclusions 


Demonstrated  by  EGRETT  measurements  that 
clear  air  and  refractive  turbulence  should  be 
sought  in  the  strong  shear  layers  found  above 
and  below  the  jet  core  of  the  winter  half¬ 
season  sub-tropical  jet  stream  in  both 
hemispheres  when  the  local  anisotropy 
represented  by  Cw^  fC^  ,the  ratio  of  vertical  to 
longitudinal  turbulence  is  weak  to  moderate: 
less  than  1.33  but  greater  than  0.6.  Future 
analyses  of  the  correlation  of  fluctuating 
velocity  and  temperature  with  fluctuating  static 
pressure  gradients  should  confirm  its  role  in 
controlling  anisotropy. 
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Temperature  Structure  Parameter 


Figure  1.  Vertical  Profiles  of  Temperature 
Structure  Parameter. 


Turbulence  Parameters 

Australia  -  6  Aug  99 


Figure  2.  Vertical  profiles  of  temperature 
variance,  turbulent  velocity  components,  and 
the  vertical  heat  flux  (multiplied  by  -1). 
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Figures.  Vertical  Profiles  of  Three 
Components  of  Reynolds  Stress. 


Table  1.  Structure  function  parameters  and  ratio 
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Table  2  Summary  of  turbulence  measurements  for  6  August  1999. 
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